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We report results of numerical and analytical studies of the spontaneous symmetry breaking in 
solitons, both two- and one-dimensional, which are trapped in H-shaped potential profiles, built 
of two parallel potential troughs linked by a narrow rung in the transverse direction. This system 
can be implemented in self-attractive Bose- Einstein condensates (BECs), as well as in a nonlinear 
bulk optical waveguide. We demonstrate that the introduction of the transverse link changes the 
character of the symmetry-breaking bifurcation (SBB) in the system from subcritical to supercritical 
(in terms of the corresponding phase transition, it is a change between the first and second kinds). 
A noteworthy feature of the SBB in this setting is a non-monotonous dependence of the soliton's 
norm at the bifurcation point on the strength of the transverse link. In the full 2D system, the 
results are obtained in a numerical form. An exact analytical solution is found for the bifurcation in 
the ID version of the model, with the transverse rung modeled by the local linear coupling between 
the parallel troughs with the 5-functional longitudinal profile. Replacing the 5-function by its finite- 
width Gaussian counterpart, similar results are obtained by means of the variational approximation 
(VA). The VA is also applied to the ID system with a mixed linear and nonlinear transverse localized 
coupling. Comparison of the results produced by the different varieties of the system clearly reveals 
basic features of the symmetry-breaking transition in it. 

PACS numbers: 05.45.Yv, 03.75.Lm, 42.65.Tg 

I. INTRODUCTION 

Symmetric double- well potentials is one of fundamental settings studied in quantum mechanics [ij and in the theory 
of optical guided- wave propagation, which obeys the Schrodinger equation similar to that known in quantum mechanics 
Counterparts of the double- well potential in optics are represented by directional couplers Jsl, l^j, including various 
dual-core waveguides created in photonic-crystal matrices [6] . It is commonly known that the ground state produced by 
the linear Schrodinger equation keeps the symmetry of the underlying double- well potential On the other hand, 
the introduction of the self-attractive nonlinearity, which transforms the linear equation into the Gross-Pitaevskii 
equation (GPE) for a Bose-Einstein condensate (BEG) of interacting atoms, loaded into the double-well potential 
Q, or the nonlinear Schrodinger equation (NLSE) modeling nonlinear dual-core waveguides in optics leads to the 
ubiquitous effect of the spontaneous symmetry breaking. As a result, the symmetric ground state is replaced, via the 
symmetry-breaking bifurcation (SBB; in fact, it is a variety of phase transitions), by an asymmetric state providing 
for a minimum of the system's energy, when the strength of the nonlinearity exceeds a critical value. This effect was 
originally discovered in a discrete model of self-trapping [8J , and later studied in many settings @ . Manifestations of 
the spontaneous symmetry breaking were also studied in detail in BEG models [l^ and demonstrated experimentally 



in a self- repulsive BEG [11 1 



In nonlinear optics, the SBB was analyzed in Ref. Q for continuous- wave (spatially uniform) states in the model of 
dual-core waveguides. For self-trapped modes in dual-core optical systems, i.e., solitons, the bifurcation was studied in 
Refs. f3,[5|. A specific example is the SBB for gap solitons in dual-core fiber Bragg gratings [l2j. Later, manifestations 
of the SBB for matter-wave solitons held in dual-trough potential traps (including solitons of the gap type, supported 
by a periodic optical-lattice potential) were explored too [l3|-[l6|. The difference of the dual-trough configuration 
from the usual double well is the presence of an additional free direction, transverse to that in which the double-well 
potential acts, see Fig. [T] below. 

A characteristic property of solitons trapped in the symmetric dual-trough potential is that, under the action of the 
self-focusing cubic nonlinearity, they exhibit the SBB of the subcritical type. In that case, branches of the asymmetric 
states emerge as unstable ones, going backward from the SBB point and getting stable after switching their direction 
forward at turning points ^V^. Thus, the system features the bistability before the bifurcation point, and this mode 
of the spontaneous symmetry change may be understood as the phase transition of the first kind. On the other hand, 
the addition of a sufficiently strong periodic potential acting along the troughs changes the character of the SBB in 
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the self- attractive medium from subcritical to supercritical. In the latter case, the asymmetric branches emerge as 
stable ones, immediately going in the forward direction [l^[l^, which may also be realized as the phase transition of 
the second kind. The above-mentioned SBB for gap solitons in the dual-core fiber Bragg grating is of the supercritical 
type too (l^ . 

Another realization of effective double- well potentials is provided by settings based on double-peak spatial modula- 
tions of the local nonlinearity coefficient, which may be implemented in optics and BEC alike, see recent review [18] of 
the topic of solitons in nonlinear potentials. The simplest form of this setting has the nonlinearity concentrated at two 
points, in the form of a symmetric pair of delta- functions or narrow Gaussians [l9j . as well as a two-dimensional (2D) 
counterpart of the system, based on the set of two parallel stripes 21 1 , or two circles [l^ carrying the self-attractive 
nonlinearity. The SBB of solitons in the double-well nonlinear potentials also features the symmetry breaking of the 
subcritical type [H, [H, [U . 

Further, the spontaneous symmetry breaking was analyzed for ID and 2D solitons in dual-core discrete systems, 
with the uniform coupling between two parallel chains (23j . or with the coupling concentrated at a single site [2^ . 
In the former and latter cases, the SBB is subcritical and supercritical, respectively. Another implementation of the 
SBB in discrete settings was recently reported for a pair of nonlinear sites embedded into or side-coupled to a linear 
host lattice [25| . 

The objective of the present work is to consider the SBB of self- attractive localized wave fields trapped in H-shaped 
potential landscapes, i.e., two parallel troughs linked by a transverse rung, as shown in Fig. [l]below. This configuration 
can be implemented in effectively two-dimensional BEC, using a set of attractive (red-detuned) laser sheets and/or 
blue-detuned repelling sheet pairs [13], or in BEC layers isolated by a strong standing optical wave [11], that can 
also be combined with magnetic trapping fields [1^. It is also possible to use the techniques allowing one to "paint" 
complex potential landscapes (in fact, even more complex than the H-shaped ones that we aim to consider) by rapidly 
moving laser beams [SOj . or induce time-averaged adiabatic landscapes created by means of variable magnetic fields 
[3l| . Essentially the same effective potentials can be created in nonlinear optics, using properly patterned photonic- 
crystal media Q , or a transverse trapping structure permanently written in bulk silica [32| . We aim to consider both 
the full 2D model with the transverse H-shaped potential (similar to the 2D model with the dual-trough potentials 
considered in Rcfs. [13, [3]) and its simplified ID counterpart (cf. the ID version of the dual trough introduced in 
Rcf. [13]). 

The transverse link (rung) added to the dual-trough potential can be used to control the dynamical properties of 
the system and, eventually, alter the character of the spontaneous symmetry breaking in the system. In particular, 
the recent results reported for the discrete system [13] suggest a possibility to change the type of the SBB from sub- 
to supercritical (i.e., the kind of the respective phase transitions from first to second) by gradually increasing the 
strength of the transverse link. 

The paper is structured as follows. The 2D and ID models are introduced in Section II, and the full 2D system is 
considered in Section HI by means of numerical methods. A set of bifurcation diagrams indeed demonstrates a switch 
from the sub- to supercritical SBB in the 2D setting. While one may expect that the strengthening of the linear 
coupling between the parallel troughs should lead to an increase of the critical value of the nonlinearity strength at 
which the SBB happens, we observe that, with the introduction of the transverse link, the critical nonlinearity strength 
at first decreases, due to the change of the character of the bifurcation, and only later starts to grow. In Section IV, we 
deal with the ID versions of the system. In that case, the transverse link reduces to a localized linear coupling between 
two one-dimensional GPEs/NLSEs. In that context, we consider different longitudinal profiles of the coupling. First, 
approximating it by a i5-function of the longitudinal coordinate, we report exact analytical solutions for the solitons 
of all the types — symmetric, antisymmetric, and asymmetric. Accordingly, an exact comprehensive solution for the 
SBB is available. Next, we consider the coupling localized in a finite interval, with the Gaussian profile. In that case, 
we develop a variational approximation (VA), and verify its predictions by comparison to numerical findings. Finally, 
we consider the ID system which combines the linear and nonlinear localized couplings between the one-dimensional 
GPEs/NLSEs, both with the Gaussian profile. For that purpose, the VA is developed too. In all the versions of the 
ID system, the SBB is always found to be of the supercritical type, unlike the full 2D model that reveals the transition 
to the supercritical type from the subcritical one. The paper is concluded by Section V. 
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FIG. 1: The two-dimensional H-shaped trapping potential. 



II. THE MODEL 

A. The two-dimensional setting 

We introduce the model in terms of the two-dimensional GPE for the mean-field wave function of the self- attractive 
BEG trapped in the H-shaped potential. In physical units, the usual form of the equation is 



2m 



(1) 



where m and a, < are the atomic mass and the respective scattering length, is the confinement length in the 
transverse direction [uz — \/'hl (mf2), if the confinement is provided by the harmonic-oscillator potential, (l/2)mf2^2:^], 
and the trapping potential of depth f/o is defined as follows: 



-f/o, at \L <\x\<D ^ \L, 
\J{x,y) = \ and at |x| < \L, \y\ < \W, 

elsewhere. 



(2) 



As shown in Fig. [U is the width of the two longitudinal troughs, with separation L between them, and W is the 
width of the transverse rung. The norm of the wave function gives the total number of atoms in the condensate, 



\^\'^dxdy. 



(3) 



Using scaled variables x = x/D, y ~ y/D, t = th/{mD'^), and ^I^ = 2y 2\/27r |as| /uzD'^, we cast Eq. ([T]) into the 
dimensionless form: 

i-^t = + -^yy) + U{x, y)m - l^-p*, (4) 

where the rescaled norm is N^d = 47rA'p|as|, and the potential takes the form of 

I -ai, at \a2 < \x\ < I + \^a2, 
U{x,y) = l at \x\<\a2, \y\ < ^a^i, (5) 



2^^^ If I ^ 2 

0, elsewhere. 



with rescaled constants 



Oil = ^2 ' "2 = "3 = ■ (6) 
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Thus, the system is governed by the set of four dimensionless parameters: norm A''2d, scaled potential depth ai, the 
relative separation between the parallel channels, and the relative width of the rung, a^. In the next section, we 
apply a numerical method to search for symmetric and symmetry-broken localized modes supported by the interplay 
of the H-shaped trapping potential and self-attractive nonlinearity in Eq. Q. We will identify the SBB in this model, 
and produce the respective bifurcation diagrams. The predicted results can be readily translated back into physical 
units be undoing the above transformations. 

In the application to the optical guided-wave propagation, Eq. (j?]) plays the role of the NLSE for the evolution of 
the amplitude of the guided electromagnetic wave, with time t replaced by the propagation distance, z. In the latter 
case, the effective potential represents the modulation of the refractive index in the transverse plane, and the cubic 
term accounts for the Kerr self-focusing. 



B. The one-dimensional system 

The ID limit of the model corresponds to the case of W" <C D <C i, i.e., as <C 1 <C a2, in terms of the relative 
parameters defined in Eq. Then, approximating the wave function in each relatively narrow trough by ID wave 
functions ■01^2 (y, t), the tails of the wave functions channeled by the narrow transverse rung in the transverse direction 
take the form of 



(V'l,2)tail ~ "^1,2 {y, t) CXp 

[which implies ai < tt^/ (2a§); otherwise (if the rung is very deep), the effective coupling between the parallel troughs 
will be stronger than given below by expression Q]. Then, using the general formalism elaborated for the analysis 
of the interaction between far separated 2D localized modes , it is easy to calculate the effective Hamiltonian of 
the interaction between the parallel troughs, mediated by the channeled tails ([7]), and thus approximate the full 2D 
model ([1]), ([5]) by the system of coupled ID equations with the local coupling: 

idtiJi = -^dyifji - lipifipi - K(5(?/)-02, 

(8) 

idti'2 = -2^y^2 - |V'2pV'2 - 

with the coupling coefficient identified by equating the above-mentioned 2D interaction Hamiltonian to its counterpart 
corresponding to the ID system ([5]): 

K = y^TT^ - 2aia| exp — -^/tt^ - 2aiQ;§^ . (9) 

In fact, by an additional rescaling of Eqs. ([5]) it is possible to fix k = 1, which is adopted below, in Section IV. To 
understand the genericity of the results produced by the ID system with the (5-functional coupling profile, we will also 
consider the ID system with the ^-function replaced by the Gaussian profile with a finite width, j/Oi 

5{y) exp i-yVvl) , (10) 

keeping the coefficient in front of the Gaussian to be 1. 



III. SYMMETRIC AND ASYMMETRIC TWO-DIMENSIONAL SOLITONS 

Stationary soliton solutions to Eqs. ([5]) were found by means of the imaginary-time integration method js^ . 
The stability of the so generated solitons was then tested by direct simulations of the perturbed evolution in real time. 
Typical examples of stable symmetric and asymmetric 2D solitons are shown in Fig. [2] 

Varying the set of control dimensionless parameters, (ai,a2,Q!3), we identified the critical value, iVcr, of the scaled 
norm A'2d, at which the SBB occurs and asymmetric states appear. The natural measure for the asymmetry of such 
states is defined as 



^2D 



(11) 



FIG. 2: (Color online) Examples of stable 2D symmetric (left) and asymmetric (right) solitons. The respective values of the 
scaled norm are N2y> = 4.05 and A'2d = 4.5, respectively. The other parameters are a\ = a2 = a-i = 1. 




FIG. 3: (Color online) Bifurcation diagrams showing the asymmetry of the 2D solitons as a function of the total norm. The 
parameters are ai = 1,0:2 = 4 in (a), and ai = 02 = 1 in (b) (fixed depth of the parallel troughs and the distance between 
them), while the scaled width of the transverse rung is varying: 0:3 = 0.5 (orange curves in each panel), az — 0.25 (blue curves) 
and as = (black curves). In the case of 03 = (no transverse linkage through the rung), the bifurcation is subcritical, 
in agreement with Ref. T4] (only stable portions of the subcritical diagram are displayed in this case). In other cases, the 
character of the bifurcation is clearly supercritical. 



As mentioned in the Introduction, the 2D system without the transverse rung, which is tantamount to the present 
2D model with W = = as well as its ID counterpart [l^, based on the system of ID equations uniformly 
coupled by linear terms, feature the SBB of the subcritical type. The difference of the present model is that, with the 
increase of 03 from zero, the character of the bifurcation quickly switches to supercritical (in other words, the kind 
of the respective phase transition switches from first to second), as shown in detail by means of the set of bifurcation 
diagrams in Fig. [31 where control parameters ai are a2 are fixed, while is varied. In terms of the underlying 
system, this means the gradual increase of the width of the transverse rung in Fig. [1] As verified by direct simulations 
(not shown here in detail), all the solution branches displayed in Fig. [3] are dynamically stable. 

The transition to the supercritical bifurcation is a consequence of the enhancement of the local transverse coupling 
between the troughs through the rung. A qualitatively similar change of the bifurcation was recently observed in dual 
discrete chains, with the transition from the uniform transverse linear coupling [l^ to that at a single site (2^ . 

The bifurcation picture is further characterized, in Figs. S] and [SJ respectively, by dependences of the value of the 
norm at the bifurcation point, A^cr, on the relative width of the transverse rung, 0:3, and the relative distance between 
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FIG. 4: (Color online) The total norm of the wave field at the point of the symmetry-breaking bifurcation, in the 2D model, 
versus the relative width of the transverse rung, 03. The other parameters are ai = 1, Q2 = 4 (a) and oi = Q2 = 1 (b). 




FIG. 5: (Color online) The norm at the point of the symmetry-breaking bifurcation in the 2D model versus the relative 
separation between the parallel troughs, 02. The other parameters are qi = 1,03 = 0.25. 



the troughs, a-i. Because, as said above, the increase of aa (making the rung wider) implies the strengthening of 
the linear coupling between the parallel troughs, one should expect the growth of N^x with 03. This is, generally, 
observed in Fig. |4l but after an initial decrease. This surprising feature is explained by the above-mentioned change 
of the character of the bifurcation, from subcritical to supercritical. On the other hand, the monotonous decrease of 
with the increase of 012 is a natural consequence of the weakening strength of the coupling between the toughs. 
The coupling between the troughs is, obviously, sensitive to the relative width of the separating barrier, with respect 
to the troughs, which is measured by see Fig. [T]and Eq. ([SJ. The fact that the bifurcation diagrams, and the 
respective critical values of the norm, are quite similar in panels (a) and (b) in Figs. [3]and|31 which pertain to a2 = 4 
and a2 = 1, respectively, demonstrates that the type of the symmetry breaking reported in this section comprises a 
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broad parametric area. 

Finally, the simulations demonstrate that the 2D solitons suffer the collapse, in the present model, exactly when it 
is expected, i.e., when the total scaled norm exceeds the well-known threshold value, A'^thr ~ 5.85 [35] . 



IV. THE ONE-DIMENSIONAL SYSTEM: EXACT, VARIATIONAL, AND NUMERICAL SOLUTIONS 
A. Exact solutions for the linear coupling with the delta-functional profile 

1. General analysis 

A remarkable feature of the ID system based on Eq. (|8]) with the transverse coupling accounted for by the S- 
functions is a possibility to find exact stationary solutions for the solitons of all the types, symmetric, antisymmetric, 
and asymmetric (note that exact solutions are not available in the discrete counterpart of the system considered in 
Ref. [24|). Stationary solutions to Eq. with a given (negative) chemical potential, fi (in the application to optics, 
— ^ is the propagation constant), are looked for as 

{^1 (y, t) , V2 (y, t)} = e-^^* My), viy)} , (12) 

where real functions u{y) and v(y) obey the following equations: 

1 tPu 3 

(13) 

[recall k = 1 is fixed in Eq. ([S]) by means of rescaling]. As follows from Eqs. ([13]), at y = the solutions must satisfy 
the following boundary conditions: 

A {u') = -2v, A {v') = -2m, (14) 
where A stands for the jump of the derivative at y — Q. Solutions to Eqs. and (fHl) are looked for as 



u = r] sech {ij {\y\ + a)) , w = 577 sech (77 (|?;| +6)), 77 = V^2/x, (15) 

where s — +1 and s = —1 correspond to the symmetric and antisymmetric states, respectively (or their asymmetric 
deformations), while shifts a > and b > are determined by equations following from the substitution of ansatz 
([T5|) into Eqs. dH]): 



V 



sinh {rja) s 
cosh^ (77a) cosh (776) ' 

(16) 

sinh (776) s 
'cosh^ (7/6) cosh (77a)' 



Remind 77 = ^/—2fi is treated here as an arbitrary constant parameterizing the family of solutions. 



2. Symmetric and antisymmetric solutions 



Symmetric modes correspond to s = +1 and a = 6, in which case Eq. ()16p yields 

exp (277asymm) = (77 + 1) / (77 - 1) . (17) 

As seen from here, the symmetric solution exists for 77 > 1, i.e., ^ < —1/2, with the total norm which can be readily 
calculated: 

r+co 

/ [u\y)+v\y)]dy = Nu + N^o^ 4 {7^-1). (18) 
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Antisymmetric states, with s = —1 and a = b < 0, are also possible, with \a\ given by the same expression ([TTl) as 
in the symmetric case. However, the antisymmetric states are expected to be completely unstable, as they correspond 
to a maximum, rather than minimum, of the respective interaction Hamiltonian, therefore they are not be considered 
below (the instability of the antisymmetric states is corroborated by numerical tests, see, e.g.. Fig. |5] below). 

3. Asymmetric solitons 

For the most interesting asymmetric solutions with a ^ b and s = +1, Eq. (|16p can be transformed into the 
following system: 

tanh^ (rya) + tanh^ (776) = 1 - 77"^ (19) 
tanh (rya) • tanh (776) = 77"^, (20) 



an explicit solution to which is 



tanh(77a) = ^(^Vl + 77 ^-^Jl-'Sri 
tanh(7;6) = i (_^\/TT^ + \/T~~3rf^^ . 



(21) 



As seen from these expressions, with the increase of 77, i.e., with the growth of the norm of the symmetric solutions 
[see Eq. (|18p ]. the asymmetric modes emerge and exist at 

77 > 7/0 = \/3, (22) 
with the following values of the norm in the two troughs: 



Nu = 277 - ^rf + 1 + - 3, 

(23) 

= 277 - ^rf + 1-^7]^ ^3. 
Accordingly, the total norm of the asymmetric solution is 

N = N^ + N, = 2 (27; - ^rf + 1) , (24) 

or, inversely, constant 77 may be expressed in terms of the total norm, which is then treated as the intrinsic parameter 
of the family of the asymmetric solitons: 




(25) 



The corresponding asymmetry measure is 



(26) 

Substituting here 77 from Eq. (j25p . we obtain an eventual expression for the asymmetry parameter as a function of 
the total norm: 



^ (iV2 - if - 12 (27V - + ^^^^ 



2 (A^2 - 4) - ^J (7V2 -4)2 + 4 (2iV - VN^ + 12)^ 

which provides a full analytical description of the SBB in the ID system. The norm of the asymmetric solitons 
assumes values N > Nq = 4 [VS— l) ~ 2.93, with Nq corresponding to the symmetry-breaking point ([2^ . 

The bifurcation diagram predicted by Eq. (|27l) . i.e., as a function of N, is displayed in Fig. [6] Obviously, the 
bifurcation revealed by the exact solution is supercritical, being quite similar to the supercritical SBB in the 2D 
model, which is displayed above in Fig. [3] Direct simulations (with the ideal (5-functions replaced by their regularized 
versions, see below) confirm that, as expected, all the solution branches shown in Fig. [5] are stable. 
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FIG. 6: (Color online) The bifurcation diagram produced by the exact solution of the ID system with the 5-functionaI profile 
of the linear coupling. 



B. The coupling with the Gaussian profile 



To check how generic the exact sohition found with the (5-functional couphng profile is, it is natural to compare the 
results to those generated by the Gaussian profile (ITU|) . Because exact solutions are not available in this case, we will 
tackle the problem by means of the VA (variational approximation), which is an effective tool for the analysis of a 
broad class of systems similar to the present one [SG^. 

To apply the VA, we use the expression for the energy of the ID system ([5]), Hr 



E = 



+ 00 



dy 



1 ^ 2 1 ^ 



i=i 



- exp 



) (VlV2*+V'l*^2) 

^0 



and introduce the following ansatz (hereafter we use ± instead of subscripts 1 and 2): 

' ,,2 



N{l±v) 
V 2W^V^ 



exp 



where W is the width of the soliton, and v is the asymmetry parameter defined as 



1 

N 



-\-oo 



dy (|V^ip-|^2p) 



(28) 



(29) 



(30) 



cf. Eq. ([Tl]) . Substituting ansatz (f29)) into expression ([SS)) yields the energy as a function of variational parameters, 
W and i^: 



E = 



^/y[+W WV2^ 

The variational equations following from expression (|31l) . dE/dW = dE/dv ~ 0, take the following form: 

1 2?/o\/l - N{1 + v^) 



(31) 



(2/0 + W^f/^ 



2/0 



N 



- 0, 



(32) 



^{yl + W^){l~v^) 2^/2nW 



10 




FIG. 7: (Color online) The bifurcation diagram produced by the ID system with the Gaussian profile of the coupling, Eqs. 
(HJ and (llOf) . in the framework of the variational approximation. The width of the Gaussian profile is j/o = 0.5 (squares), 0.25 
(stars) and 0.125 (balls), respectively. 



These equations were solved in a numerical form, to find parameters W and v of stationary solutions as functions of 
the free parameters, yo and N . The stability of the solutions was verified, in the framework of the VA, by checking if 
they correspond to a local minimum of energy (1311) . 

The so generated bifurcation diagrams are displayed in the left panel of Fig. [7] for several fixed values of width ijq 
of the Gaussian profile of the coupling. It is clearly seen that the SBB is supercritical in this case too, being quite 
similar to that produced by the exact solution for (5-functional profile, cf. Fig. [SJ and to the supercritical diagrams 
found in the 2D model, cf. Fig. [3] The respective dependence of the norm at the bifurcation point on the width of 
the Gaussian profile, iVcr(yo)- In fact, the latter plot is a ID counterpart of the one which, in the 2D model, shows 
the dependence of Ncr on the relative width of the transverse rung, as, cf. Fig. S) 

We have also constructed stationary states of the ID system with the Gaussian-shaped coupling as numerical solu- 
tions of Eqs. ([U and (fTO|) . Two different algorithms were used for this purpose, the imaginary-time integration and the 
Newton-Kantorovich iteration method, both yielding identical results. Thus, three types of stationary configurations 
were found — symmetric, antisymmetric, and asymmetric ones. Their stability was tested by direct simulations of the 
evolution in real time. It has been confirmed that the respective SBB is indeed supercritical, and the variational re- 
sults are very close to the numerically exact ones. It was also found that antisymmetric modes (which do not undergo 
the bifurcation) are unstable (as mentioned above), see an example in Fig. |8l The two-peak profile of the stationary 
mode in this figure is a consequence of its antisymmetry (the Gaussian profile of the coupling is effectively repulsive, 
in this case, which also explains the instability of the mode). It is seen that almost the entire norm is spontaneously 
transferred into a single trough, where the two peaks split into separating single-component solitons. 

We have also developed a similar analysis, using both the VA and numerical solutions, for the ID model with a 
rectangular profile of the transverse coupling. The results (not shown here) are very similar to those generated by the 
Gaussian profile. 



V. THE COMBINED LINEAR-NONLINEAR COUPLING WITH THE GAUSSIAN PROFILE 

It may also be relevant to take into regard a nonlinear correction to the linear coupling between the two troughs, 
in the framework of the ID system (note that the full 2D model automatically takes into regard both linear and 
nonlinear effects of the coupling). To this end, we adopt the following modification of the ID system (|8]), with the 
Gaussian coupling profile (flU)) : 

exp (-2/V2/0) (o-|V'2pV'i + V'2) , 

(33) 

exp (-yVyo) {'^li'i + V'l) , 
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FIG. 8: (Color online) The shape of a typical antisymmetric solution and its evolution in the ID system ([8]) with the Gaussian 
profile (|10p of the linear coupling. In the top panel, the pink, red and (dash) blue curves display tpi (y) , xf}2 (y) in the stationary 
solution, and Gaussian profile exp(— y^/j/o), respectively. The two lower panels show the evolution of {y,t)\ and \'ip2{y,t)\, 
respectively. The parameters are: total norm A'^ = 11.4, chemical potential /i = —1.03, and the width of the Gaussian profile 
2/0 = 0.5. 
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FIG. 9: (Color online) The same as in Figs. [7] , but for the ID model with the Gaussian profile of the mixed linear-nonlinear 
coupling, based on Eq. (|33|l and the variational ansatz (|29|) . The width of the Gaussian profile is yo — 0.25 (squares), 0.125 
(stars) and 0.05 (balls), respectively. The relative strength of the nonlinear and linear couplings is ctq ~ 1- 



where a is the relative strength of the nonlinear coupling. In this case, the energy functional is 



E 



■ exp 



/ + OC I 1 ^ 1 ^ 



2/5 



(34) 



To apply the VA method, we adopt the same Gaussian ansatz (P^)) which was used above. The substitution of the 
ansatz into functional ([34]) yields 



N 



Nyo^ 



l-iy2 



4VK2 "''\]W^ 



the corresponding variational equations, dE/dW = dE/dv = 0, being 



(35) 



(W2+y2)3/2 

^i:{W^ + 2yl) \ + 2yl 



1 



= 0, 



(36) 



2yo 



N 



V(VK2+2/2)(i_ j,2) \ + 2yl 



The bifurcation diagram produced by the numerical solution of Eqs. (1551) . and the respective dependence of the total 
norm at the SBB point on y^ are displayed (for g — V) in Fig. [5] 

The comparison of the dependences Ncr{yo), obtained in the three versions of thelD model with the finite width 
of the coupling profile (Gaussian, rectangular — which is not displayed here in detail — and mixed linear-nonlinear) 
suggests that the inclusion of the nonlinear coupling into Eq. ([55]) makes this characteristic essentially more similar 
to its counterpart obtained in the full 2D model, cf. Fig. |4l While in Fig. [71 and in its counterpart obtained in 
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the model with the rectangular profile, the curves Nci{yo) are convex, they are concave in both Figs. |4]and|9l This 
conclusion is naturally explained by the above-mentioned fact that the full 2D model automatically takes into account 
both the linear and nonlinear coupling between the parallel troughs. 

VI. CONCLUSION 

The objective of this work is to extend the study of the spontaneous symmetry breaking in ID and 2D solitons 
trapped in the systems of two parallel tunnel-coupled potential troughs. Unlike recently studied models, here we 
introduce the H-shaped potential, and focus on the SBB (symmetry-breaking bifurcation) for solitons in this system. 
The analysis of both the 2D and ID versions of the model demonstrates that the transverse link (the rung of the 
H-shaped potential profile) transforms the bifurcation from subcritical into the supercritical one. In other words, the 
rung controls the switching of the corresponding symmetry-breaking phase transition from the first into second kind. 
A nontrivial manifestation of the change of the SBB type is the non-monotonous dependence of the critical value of 
the soliton's norm, at the bifurcation point, on the strength of the transverse link: prior to the expected growth, it 
demonstrates a region of the decrease. In the full 2D model, the results were obtained in the numerical form. On 
the other hand, the ID version with the (S-functional profile of the local linear coupling between the troughs admits 
the exact analytical solutions for the solitons of all the types — symmetric, antisymmetric, and asymmetric — and, 
accordingly, the bifurcation diagram was obtained in the exact analytical form, which is a rare possibility. In other 
variant of the ID model, with the Gaussian profile of the local transverse linear coupling, the results were obtained 
by means of the VA (variational approximation). The VA was also used to solve the bifurcation problem in the 
system combining the linear and nonlinear localized transverse coupling. The set of the results reported in the paper 
provides for a comprehensive description of the symmetry breaking in the H-shaped system. In particular, adding the 
nonlinear coupling to the ID model makes the characteristics of the symmetry breaking closer to those found in the 
full 2D model, which automatically incorporates linear and nonlinear coupling effects. 

The settings considered in this work can be realized in the self- attractive BEG, and in self- focusing bulk optical 
media, with the appropriate transverse profile of the refractive index. It is relevant to stress that the description of 
the BEG based on the coupled GPEs is entirely based on the mean-field approximation. While it is well known that 
this approximation provides for an exceptionally accurate description of virtually all matter- wave patterns observed in 
experiments, quantum fluctuations and other beyond-mean-field effects being detectable only under special conditions 
0, one may expect that fiuctuations can be amplified in a vicinity of the phase transition (i.e., of the SBB), which 
suggests a question, whether the GPEs are applicable in this case. A full analysis of this issue should be a subject for 
a separate extended work; nevertheless, a relevant analogy is suggested by the analysis of fluctuations in the vicinity 
of the SBB, both in the continuous-wave and soliton settings, which was performed in the framework of the model of 
the directional coupler in nonlinear optics, based on the coupled ID NLSEs, in Ref. Q. The conclusion was that, on 
the contrary to the anticipations, the fluctuational effects remain virtually negligible in that case, their amplification 
in the vicinity of the respective SBBs producing no conspicuous changes against the mean-field description. Based on 
this analogy, we expect that the mean-field description remains valid in the present case too. 

One possibility for further development of the analysis is to perform it in a two-component system. Another 
interesting extension may be carried out for a configuration with two parallel transverse rungs, in which case one may 
expect a double spontaneous symmetry breaking: between the parallel potential troughs, and, independently, between 
vicinities of the two rungs. In particular, a challenging problem would be to find an exact solution in the case when 
the pair of the rungs is represented by a symmetric set of two (5-functions. 
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